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Abstract 

Scalar fields at finite temperature are considered in four dimensional ultrastatic curved spacetime. 
One loop nonlocal effective action at finite temperature is found up to the second order in curvature 
expansion. This action is explicitly infrared finite. In the high temperature expansion of free energy, 
essentially nonlocal terms linear in temperature are derived. 
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1 Introduction 



Effective action is a powerful tool in quantum field theory. It contains information about all 
Green functions, vacuum polarization effects, etc. In general case it is a non-local functional. 
To study thermodynamic properties of a system one has to know the effective action at non- 
zero temperature. In this letter we calculate the one loop nonlocal effective action and 
free energy for quantum scalar fields in a curved spacetime at finite temperature T. At 
first we introduce general notations for both zero and finite temperature theories on some 
Z)- dimensional manifold. We consider fields ip that satisfy the equation 

F(V)¥^ = 0, F(V) = n + P(x)-^i?(x). (1) 

o 

Here the Laplace-Beltrami operator □ = fi'^'^V^Vi, is constructed in terms of a covariant 
derivative V^. The metric g^^ is characterized by its Riemann and Ricci curvatures R^.^vp = 
duT'^p — • • • , Rai3 = R^afii3- The term with scalar curvature R/6 in the operator -F(V) is 
explicitly singled out for convenience only. The potential P{x) has an arbitrary structure 
and may include a curvature R{x) or an interaction potential. This means that models 
satisfying eq. (|l|) are quite general and describe non-conformal and interacting scalar fields. 

The one loop effective action can be expressed in terms of eigen values of the operator 
-F(V) and formally reads W = ^TrlnF. It can be written || in the form of an integral 
over the proper time s, 

1 d s 

W = -- -TtK{s). (2) 

2 Jo s 

In coordinate representation Tr denotes the functional trace Tr K{s) = J d^xK{s\x,x) of 
the heat kernel, which is defined as a solution of the equations 

^K{s\x, y) = F{W)K{s\x, y), K{0\x, y) = 5{x - y). (3) 

ds 

Knowledge of the heat kernel with separated points allows one to construct the covariant 
diagrammatic expansion of the effective action to all loop orders For our case of one 

loop computations, it is sufficient to find the trace of the heat kernel. 

Thermodynamics of a system is well defined if background fields are static, so we restrict 
ourselves to considering static spacetimes. It is well known p, ||, Q that a quantum field 
theory defined in a Euclidean spacetime, which is periodic (for bosons) in imaginary time 
r with a period /?, amounts to a thermal field theory with temperature T = 1/13. In this 
letter we consider Euclidean manifolds with the topology of a cylinder xR^. The periodic 
(thermal) heat kernel can be expressed as an infinite sum of zero temperature (vacuum) 
heat kernels |^ 

oo 

K'^{s\t,x;t',x') = K{s\t,x;t' + I3n,x') (4) 

n=— oo 

This image sum is equivalent to summation over Matsubara frequencies in a momentum 
space representation in thermal field theory. 
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We restrict now a class of spacetimes to the one possessing ultrastatic (optical) metrics 

ds^ = g^^dx^'dx" = dr^ + gij{x)dx''dx^ (5) 

where /x, = (0,1,2,3) and i,j = (1,2,3), and the metric tensor gij{x) depends on the 
spatial coordinates x only. To take into account temperature effects in curved spacetime 
explicitly one has to factorize the heat kernel into temporal and spatial K^^^ (s) parts, 

Kis\T,x;T',x') = exp(- "j'^' )K^'\s\x; x'), (6) 

what is possible to do in ultrastatic spacetimes. Here and below we denote objects in 
ultrastatic spacetimes by overbars. Then, the trace of (^ takes an elegant form |^, |^ in 
terms of Riemann theta function, which is defined in a usual way 6*3(0, h) = Y^n='^oo e^"'*'6'^ , 

TrKP{s) = es{0,e-^)j^^ J d'x K^'\s\x,x). (7) 

2 Nonlocal effective action and free energy at finite 
temperature 

Various methods are applicable for calculation of Tr^. A traditional tool of quantum field 
theory in curved spacetime, the Schwinger-DeWitt series |jl|, corresponding to the heat 
kernel expansion at small proper times s, is useful only for obtaining high temperature 
expansion of effective action ||^. We resort to a method of covariant perturbation theory 
P, [ly, [Tl|, , since local representations for the heat kernel are not adequate to our problem. 
This technique, as proposed by Vilkovisky in ref. |T^, is a regular procedure for covariant 
expansion of K{s) in powers of curvatures (field strengths). This is essentially nonlocal rep- 
resentation, and it takes into account an infinite number of derivatives acting on curvatures. 



At small s it reproduces the local Schwinger-DeWitt expansion [|T^, . 

Besides, the Schwinger-DeWitt expansion, which is very effective in dealing with ultravi- 
olet divergencies of the effective action [||, fails in the infrared (large s) regime for massless 
theories. In contrast, the one loop effective action and all Green functions computed in the 
framework of covariant perturbation theory are infrared finite if the dimension of a manifold 
under consideration is larger than two |P, 0, |14|. In this letter we show that the one loop 



nonlocal effective action for arbitrary scalar fields calculated on a four dimensional mani- 
fold remains infrared finite at finite temperature. This fact is very important, for infrared 
divergencies are one of major obstacles in thermal field theory 

It is necessary now to give a brief description of covariant perturbation theory at zero 
temperature in arbitrary dimensions D. Let us first display a structure of the nonlocal trace 
of the heat kernel [|T0[ , 

TtK{s) = (4 J)D/2 / d""^^'^' {l + sP + [R,M-sa)R'^'' 

+ Rf^{~sn)R + pf,{-sa)R + pf,{-sn)p] } + 0[^% (8) 
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Form factors fi are analytic functions of the dimensionless argument sD. They act on 
tensor invariants constructed of the set of field strengths R"^, P (called here curvatures) 
characterizing background. We use for these curvatures the collective notation 3?. First two 
terms of the sum (Rf) are purely local and coincide with first two coefficients of the short time 



expansion [T^ . The Euclidean spacetime is supposed to be asymptotically flat and have the 
topology R^. Therefore, gravitational curvatures and potential P should vanish at infinity 
[T0[| . The calculations in covariant perturbation theory are carried out with accuracy 0[9ft"], 
i.e., up to terms of n-th and higher power in the curvature 3ft. Thus, this curvature expansion 
is valid for highly inhomogeneous background fields, VV3ft >> 3ft^. 

All form factors in can be expressed in terms of one basic form factor 



f(-sn) = /dae°(^"°)^°. 
Jo 

Their explicit form reads [[l^ 



(9) 



/i(-^n) 

f2{-sD) 



(sn)2 

36^^ ^3 sD {say 



(10) 



Let us return now to the finite temperature case. The free energy for scalar fields in ultra- 
static spacetime is defined in terms of the Euclidean effective action at finite temperature 
by the relation 



(3 2/3 io s ^ ' 



'ID 



It is convenient to separate a vacuum mode n = from the infinite sum in (^, since 
only vacuum effective action contains ultraviolet divergencies [|17|. This is equivalent to 
subtracting the zero temperature free energy from F^. Since only renormalized effective 
action and free energy have physical meaning, we renormalize one loop vacuum effective 
action (^ with the use of zeta function regularization: 



oo 
rcn 



ld_ 
2d^ 



2e 



T(e) 



ds 



Tr/ST" 



(12) 



e=0 



where is a mass-like regularization parameter and F is the gamma function. will be 
combined with n 7^ terms at the end of our derivations. Thus, we compute 



FL - F. 



00 

ren 



1 r°°ds 



2 Jo 



(63(0, 



e 4s 



1 



1 



(47rs)i/2 



six: X] 



(13) 



The three dimensional heat kernel K^^\s) is defined as a solution of (^) with the three 
dimensional operator (|1|). In this case the operator □ is the three dimensional Laplacian 
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A, and the curvatures 3? = {P{x), Rij{x)) are defined on a three dimensional hypersurface 
r = const. 

The integral (0) with first two terms of the trace of the heat kernel (H) gives just 
numerical coefficients. This local contribution to free energy is known & [1^ and coincides 



with first two terms of high temperature expansion [|1^]. The problem reduces to computation 
of the thermal form factors 



y 



Ay 



A 



(14) 



where y = /3^/4s. 

Calculation of ([T^) is a key technical part of this letter and we explain its idea on the 
example of the basic form factor (P) . Straightforward integration over y leaves us with the 
modified Bessel function of the second kind 
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4 f; J^' da Ko {nl3^a{l-a).pK^ . 



(15) 



Another change of variables, x = 2^ a{l — a), allows us to find its form in terms of the 
exponential integrals, 



dx 



X 



^0 



,nz X, 



nz 



Ei(I^)e-"^/2 - Ei( - — )e"^/2 



(16) 



where z = A. Now we can use for the right hand side of ([16D its standard form in terms 
of elementary functions and obtain 



T 



-A) 



sin(t) 



t^ + {nzy/4' 

This is just what we need since the sum over n can be exactly evaluated. 



(17) 



^^t^ + n^zyA 2 



2tt 



zt th(27rt/z) 



1 

t2 



There are, as seen from (|TT]), two other types of basic thermal form factors with one and two 
subtractions. Derivations for them are similar but more technically involved. 
Finally, we arrive at the following integral representation. 



FL - K 



oo 
ren 



327r2 



(Tx q^'^ \ H -P 



+ 



i?,,7r(-A)i?^^+i?72^(-A)i? 



+ p4{-A)R + P^I{-A)P 
Here the thermal form factors have a form 



7f (-A) = Tdt g,{t) 
Jo 



2-K 



zt th(27rt/z) 



1 
1? 



(19) 



(20) 
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and Qi (i = 1, ...4) are the simple combinations of elementary functions 

-cos(t) sin(t)- 



92it) 

93it) 



1 /sin(t) 
2V^ + ^ 



t3 



t4 



— (-sm(t) 
48^3 ^' 



o Cos(t) ^ sin(t) 



t 



1/1 , , cos(t) 

-f-sin t + ^ 

2^3^^ t 

sin(t). 



t2 
sin(i(:) 



cos(t) sin(t) 



t3 



t4 



(21) 



We have achieved our goal of evaluating the image sum (|^), moreover the answer is 
expressed in terms of elementary functions. As can be seen from eq. ( pOD 'a golden rule' of 
field theory calculations [0, namely, to keep a proper time integration to the very end, has 
been also satisfied, as the variable t (which is dimensionless in contrast to the proper time s) 
clearly plays a role of the proper time. One can check that integrals of this type are regular, 
and two apparent singularities at the lower limit compensate each other. The price for all 
this is that an argument of the hyperbolic tangent function is the dimensionless combination 

The final result for renormalized free energy at finite temperature F^^^ is presented by a 
sum of eqs. (p!9|)-(|2T|) and renormalized free energy at zero temperature F^^. After the zeta 
regularization (|12]), the latter one takes the form, 

^ ^d^a: g^/^ {i?,,7i(-A)i?^^' + /?72(-A)i? 



327r^ 

+ P73(-A)i? + Pj^{-A)P + 0[^^ 
where form factors 7j(— A), i = 1, ...4, are 



(22) 



7i(-A) 



1 

60 



■ / A^ 46 



15 



73(-A) = - 



72(-A) = 



74(-A) 



1 



-In 



Ax 

-)+2 



2 / 



(23) 



3 High temperature expansion of free energy 

Let us emphasize that the formulae (P!^-(PT|) are valid at arbitrary finite temperature. 



Unfortunately direct evaluation of integrals ( pOD is possible only if one studies asymptotic 
behavior of the free energy, e.g., high temperature regime. It can be done relatively easy, 
because, as usual, the problem boils down to finding /5 ^ asymptotic of the thermal 
form factors (pO|)-(|2T|). One has to be careful while dealing with mutually compensating 
singularities in 7f . 
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-A) 



After some algebra, the outcome for, e.g. 
27r2 ,13^/^. 



IS, 



C(3) 
247r2 



+ 2 ln( 



C 



/5'(-A) + 



C(5) 
6407r^ 



(24) 
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where C is the Euler constant and C is the Riemann zeta function. 

Having (p!9D expanded in temperature series, we can add the vacuum free energy 
it. The result for the high temperature expansion of the renormahzed one loop free 
takes a form, 



energy 



PL 



1 



+ 



32/3 



90/^4 



24/32 



1 ^ 



12 



1152 

p 



R 



1 ^ ^ 1 ^ ^1 



—R.^R^" RR + -PP] ] 

60 1 sn 2 J J 



180 



(25) 



First and foremost, one should note that 1/P order has now form ^ i — Therefore it 



A 



is essentially nonlocal, and techniques based on local short time expansions can not generate 



it. We expect that terms of higher orders in curvatures [0, |14[ will also give nonlocal 
contribution to the linear T order. The meaning of these structures can be understood from 



an elegant spectral representation [11, 12 1 in terms of massive Green functions. 



1 



1 



2 

— — / (^Tfi 

n Jo m2 - A ' 



(26) 



A remarkable property of the expression ( pSj) is that all ln(— A) have disappeared, leaving 
logarithmic temperature dependance in the form of ln(/3/i), the combination known in flat 
space thermal field theory [1^, ^ . 

Combination of quadratic in curvatures terms at the logarithm is nothing but the trace 
of the second Schwinger-DeWitt coefficient 02, taken with Riemann curvature expressed via 
Ricci one [0, ^ |10|- Higher powers of P (see eq. (plD), which we discarded in Ff^^, are also 



quadratic in curvatures parts of higher Schwinger-DeWitt coefficients M . 

All local terms in our result (^) perfectly reproduce those of ref . |jl8[ . We also obtained 
the explicit form of all nonlocal terms of the second order in curvatures. They contain form 
factor 1 / y/—A and turn out to be proportional to 1/(3. Implicitly these nonlocal parts are 
contained in the term denoted by C(3)(0) in eq. (17) of jlSl . 



4 Conclusion 

The main result of this paper is the one loop nonlocal free energy for scalar fields at finite 
temperature, eqs. ([T9|)-(p3D. Since in the setting of our problem and differ only by 
a factor (3~^, one can speak about free energy and effective action interchangeably. One 
can derive other quantities, such as Green function, stress tensor or entropy, by varying this 
free energy correspondingly with respect to the potential, metric or temperature. Clearly, 
to proceed to physical applications one first needs some kind of spectral representation for 
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finite temperature form factors, like tlie ones derived for zero temperature effective action 

a In 



11, 12 1 or its liigli temperature limit 



It should be emphasized that the perturbation theory we used to derive our results works 
well for highly inhomogeneous background fields, i.e., VV3fJ >> 3?^. Therefore, eqs. (|19|)- 
(^), (|25|) are not apphcable to the opposite case of constant or nearly constant background. 

Our work also demonstrates, that infrared divergencies in finite temperature field theory 
are artificial and a correctly computed thermal effective action is infrared finite at all stages 
of calculations, what supports the main idea of ref. . 

Generalization of our derivations to arbitrary static spaces does not seem to pose serious 
complications, and involves well elaborated techniques of conformal transformations to ultra- 
static spacetimes [pl|, |l3]. In order to include into consideration manifolds with boundaries. 



a major revision of the perturbation theory [jTOl is needed, though it is still feasible to do. It 
is also possible to include the chemical potential Non-Abelian gauge fields can be taken 
into consideration as well, but the factorization formula will be modified by additional 
terms p2|, connected to temporal components of gauge fields. One can also easily generalize 



the results to quantum fields of other spins. 
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